Direct evidence for an intrinsic square vortex lattice in the overdoped high- Tc 

superconductor Lai 83Sro.i7Cu04+5. 
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We report here the first direct observations of a well ordered vortex lattice in the bulk of a La2-2,Sr2,Cu04+i 
(La214) single crystal (slightly overdoped, x=0.17). Our small angle neutron scattering investigation of the 
mixed phase reveals a crossover from triangular to square coordination with increasing magnetic field. The 
existence of an intrinsic square vortex lattice has never been observed in high-temperature superconductors 
(HTSC), and is indicative of the coupling of the vortex lattice to a source of anisotropy, such as those 
provided by a d-wave order parameter or the presence of stripes. 

PACS numbers: 74.60. Ge, 74.72.Dn, 61.12.Ex 



Apart from the unusual electronic and magnetic be- 
haviour of the cuprate HTSC, experiments reveal a 
tremendously rich variety of mesoscopic phenomena as- 
sociated with the flux vortices in the mixed state [Q. 
Due to their two dimensional electronic structure, the 
HTSC are highly anisotropic. The anisotropy is charac- 
terized by the ratio 7 = Aj^/A|| , where Aj_,A|| are the 
superconducting penetration depths for currents flowing 
perpendicular and parallel to the two-dimensional Cu02 
planes. In La214 the degree of anisotropy (7 ~ 20 for 
x=0.15) lies between that of the YBa2Cu30a; (Y123) and 
Bi2Sr2CaCu208+a; (Bi2212) materials [1). The cuprates 
are also extreme type-H superconductors, indicated by 
the high value of the Ginzburg-Landau parameter k — 
A/^ , where ^ is the superconducting coherence length. 
In HTSC the combination of high transition temperature 
Tc, high 7 and high k leads to exotic vortex behaviour, 
such as the phenomenon of vortex lattice (VL) melting 
101 . For a conventional (isotropic) pairing mechanism 
such anisotropic conduction properties can lead to dis- 
tortions of the vortex lattice as the applied field is tilted 
towards the Cu02 planes, but the local coordination re- 
mains sixfold 

From the magnetic point of view the single-layer (one 
Cu02 sheet per unit cell) La214 compounds appear to be 
rather different from their two-layered Y123 and Bi2212 
counterparts. La214 is characterized by the presence of 
incommensurate spin excitations located in the vicin- 
ity of the antiferromagnetic wavevector (tt, tt) of the un- 
doped parent compound and at low temperature a spin- 
gap develops in the spectrum. While such excitations 
could originate from coherence effects in the supercon- 
ducting state Q it has also been proposed that they 
could indicate the presence of dynamical stripes S] . It 



is still unclear to what extent the presence of dynami- 
cal stripes and/or antiferromagnetic fluctuations affects 
the vortex state. Recently, inelastic neutron scattering 
experiments in a magnetic field have revealed sub-gap 
excitations which have been attributed to the existence 
of antiferromagnetic excitations in the vortex core [ p^ . 
It was also shown that the opening of the spin-gap is 
correlated to the melting-line of the VL rather than 

These results indicate that a subtle interplay 
may indeed exist between the unusual microscopic and 
mesoscopic properties of the La214 compounds. 

Our experiments were performed on the small angle 
neutron scattering (SANS) instruments of both the Paul 
Scherrer Institute (PSI) Q, Switzerland, and of the In- 
stitut Laue Langevin (D22), Grenoble, France, with a 
neutron velocity selector of 10% full width at half max- 
imum, to produce a flux of neutrons with a mean wave- 
length Aq between 6 and 15 A. The incident neutrons 
were coUimated over a distance of 18 m before reaching 
the sample. A 96x96 cm^ multidetector ||l^ at a distance 
of 18 m was used to detect the scattered neutrons. 

The sample was a single crystal cylinder of length 
20 mm, diameter 6 mm with the c-axis perpendicular to 
the cylinder. One of the {1, 1} in-plane axes, which cor- 
respond to the directions of the Cu-0 bonds of the Cu02 
planes lay roughly along the cylinder axis. The crys- 
tal had composition Lai.83Sro.i7Cu04 (Tc=37 K). The 
width of the superconducting transition, measured by ac- 
susceptibility, was quite small (ATc=1.3 K, as defined by 
the 10%-90% criterion), indicating the high sample qual- 
ity. It was mounted in a cryostat in a magnetic field of 
0-1.2 T applied parallel to the incident beam. 

For the initial measurements, the c-axis was oriented 
along the field direction. In a SANS experiment a con- 
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ventional VL should give rise to Bragg reflections at 
reciprocal lattice points qhk- The intensity of a sin- 
gle (h,k) reflection is given by hk oc \Fhk\'^/qhk, where 
Fhk is a spatial Fourier component of the field profile 
in the vortex lattice. For the case of a high-K super- 
conductor such as La214 with the field perpendicular to 
the superconducting planes, in the London limit (fields 
well below the upper critical field) the form factor is 
Fhk = B/{1 + ((7/ifeA|j)^), where B is the average mag- 
netic induction. This expression may easily be modified 
to account for an anisotropic system For B > noHd 
{Hci is the lower critical field) the intensity of a first order 
reflection reduces to Iio oc dio<&o/'^| where qio = 27r/dio 
and $0 is the flux quantum h/2e. In La214, the London 
penetration depth is similar to that measured in Bi2212 
I psf . Although SANS has been used extensively to mea- 
sure VL structure in both the Bi2212 ||l6| and Y123 sys- 
tems [^-^, to date the La214 family remains virtually 
unexplored ]l7[ |. 

Figure 1 shows the difference of D22-data taken af- 
ter coohng from T=40K> to T=1.5 K at 0.1 T and 
a background taken at T=40 K. Figure 2 shows simi- 
lar PSLdata at 0.8 T with a background taken at low 
temperature in zero-field. These subtractions are a clean 
way to remove the large background signal arising from 
small angle defect scattering from the crystal. Due to 
the large mosaic of the vortex lattice, several spots could 
be observed simultaneously. We could furthermore ro- 
tate the cryomagnet and sample about the vertical and 
horizontal axes, allowing us to collect the integrated mag- 
netic intensity as we rocked the Bragg spots through the 
Ewald sphere. Figure 3a shows the tangential average of 
the neutron signal, as a function of the modulus of the 
wavevector q. As expected, the position (in reciprocal 
space) of the peak maximum changes with the field, thus 
clearly establishing the VL origin of the neutron signal. 

We now return to the lowest-field data (B=0.1 T) 
shown in Figure 1. When the c-axis lies along the field 
direction and the field is kept constant during cooling, 
we observe a large number of spots distributed around a 
ring (Fig. la, high-resolution mode: Aq = 14.5 A). Since 
the value of the wave vector obtained from the tangen- 
tial average is close to the one expected for a triangular 
lattice (see Fig. 3b), the pattern most likely consists of 
a superposition of diffraction from various domain orien- 
tations of triangular coordination. In order to confirm 
this triangular coordination, it is necessary to reduce the 
high degree of degeneracy apparent in Figure la. This 
was achieved by rotating the c-axis (about the vertical 
axis) away from the field direction by a sufhciently large 
angle 0. Further improvement of the mosaic spread of 
the VL was obtained by oscillating the field around its 
mean value while cooling. The results of these proce- 
dures are shown in Fig. lb (0 — 10°, Aq = 8 A) where a 
triangular pattern can be observed. Our measurements 



shown in Fig. 1 are very reminiscent of what has been 
seen in untwinned Y123 ||] crystals. In that system, a 
small amount of residual twin planes appears to control 
the orientation of the VL, giving rise to a high level of 
degeneracy. As the c-axis is rotated away from the field 
direction, the triangular symmetry becomes evident. 

As the field is increased from 0.1 T to 0.8 T a com- 
pletely different pattern emerges, with all of the mag- 
netic scattering concentrated in four intense spots ap- 
pearing along the {1,1} directions, forming a perfect 
square (Figure 2a). To establish the significance of this 
fourfold diffraction pattern at higher field we must first 
take account of the influence of the twin planes. While 
the alignment of such patterns has previously been at- 
tributed to the influence of an anisotropic vortex core in 
Y123 [[l8| , more recent measurements on untwinned Y123 
crystals Q have confirmed earlier suggestions ^,|l^ that 
these results in twinned crystals are actually explained 
in terms of strong alignment of the triangular VL mainly 
due to twin-plane pinning. Thus in Y123 when the field 
is roughly perpendicular to the ab planes, an intrinsic 
triangular coordination may give rise to a predominantly 
fourfold diffraction pattern. This conclusion is also in- 
ferred in twinned Y123 by rotating the field away from 
the c-axis so that the influence of the twins is severely 
reduced, resulting in the recovery of a triangular coordi- 
nation g. In Y123 for smafl tilt angles 9 < 6^ « 5° a 
fraction of vortices remains pinned to the c-axis, since the 
vortices bend in order to lie for part of their length within 
the planar defects . For larger angles this is no longer 
energetically favourable and the vortices lie along the di- 
rection of the applied fleld |^. Estimates of this critical 
angle Qc for La214 would not suggest a significantly dif- 
ferent behaviour from Y123, given the slightly increased 
values for A and 7, although the precise value is difficult 
to predict [|] and falls with increasing field. Our measure- 
ments at low field (see Fig.l) indicate that for our La214 
compound Oc is lower than 10°. Hence, for data taken at 
high field, we can be sure that O — 30° exceeds the crit- 
ical angle. We see in Figure 2b that a four-spot pattern 
is retained at this angle. Thus it is extremely unlikely 
that this pattern arises in La214 from pinning distor- 
tions due to the presence of twin boundaries, and these 
data represent the first such observations of an intrinsi- 
cally square VL in a cuprate HTSC. For completeness 
we note that for O = 30° there exists a slight distortion 
of the square pattern which is in accord with expecta- 
tions for fields tilted away from the c-axis. Due to the 
approximately uniaxial superconducting anisotropy, the 
tilting causes an increase of the penetration depth of the 
vortex currents along one direction, distorting the vor- 
tex shape and hence the vortex lattice. In the London 
approach for anisotropic uniaxial superconductors pO| a 
rhomboid is expected. The expected values for the ratio 
e = a/b between the sides (a and b) of the rhomboid and 
for the internal angle (3 are e=0.88, /? = 86.4°. A two 
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dimensional fit of the neutron data yields e = 0.90(3) 
and (i = 87.3°(1.7°), which agree well with the theoreti- 
cal values. Furthermore, the scattered intensity changes 
roughly as cos^(0) as expected from the changes in the 
penetration depth. This indicates that all of the vortices 
lie along the field direction and not along the c-axis, that 
is, they are not pinned within the twin planes. Further 
contributory evidence for a square coordination comes 
from the positions of the Bragg spots in reciprocal space. 
The relationship between the magnetic field B and the 
magnitude q of the wavevector depends on a structure- 
dependent quantity a — {2t: / q^B /<^q, where a is equal 
to a/3/2 or 1 for triangular and square lattices, respec- 
tively. At high magnetic field, the experimental values 
of a are as expected for a square lattice, while at low 
field they are consistent with that of a triangular lat- 
tice (see Figure 3b). An alternative way to quantify the 
triangular to square transition consists in monitoring the 
intensity ratio of sectors (±15°) containing the {1, 0} and 
{1, l}-type directions. As shown in Figure 3b this inten- 
sity decreases steadily from 1 at B=0.1 T to about 0.2 at 
B=0.4 T and then remains constant. 

The experimental evidence thus indicates that an in- 
trinsic square vortex lattice exists at fields 0.5 T and 
larger. We now discuss the possible origins of this four- 
fold coordination. Isotropic London or Ginzburg-Landau 
(GL) theories predict a triangular VL in high-K systems 
such as the cuprates. This is also true for modified the- 
ories which take into account the anisotropy of layered 
materials such as the HTSC. VL of square coordination 
have however been observed in non-cuprate superconduc- 
tors and may occur for a number of reasons. In 1ow-k 
materials non-local electrodynamic effects are expected 
to be more significant due to stronger interactions be- 
tween the cores of neighbouring vortices. This may give 
rise to a VL of square coordination, although at low mag- 
netic induction (large vortex separation) or close to Tc 
(reduced infiuence of non-local effects) the VL tends to a 
triangular coordination. In combination with anisotropic 
electronic structure, this explains the recent observations 
of triangular to square transitions in the nickel boron- 
carbide superconductors . A square VL can also arise 
in extended GL theories which allow for the existence of 
more than one order parameter, which may account for 
the robustness of the square VL in the candidate p-wave 
system Sr2Ru04 ||2^. A fourfold symmetry could also 
result from the anisotropic (d-wave) nature of the super- 
conducting gap p3|-p6t via the increasing importance of 
the anisotropic vortex cores at high field. The onset of 
such effects might be expected to occur at fields a factor 
of order ten lower than in Y123 (T™"^ = 93if), due to 
the smaller gap in La214 (T™'^^ = 38.5/i:). The experi- 
mental search for such a transition has been complicated 
in Y123 by the influence of twin-plane pinning discussed 
above. To date there are no reports on either twinned 
or untwinned Y123 for such an intrinsic change of coor- 



dination. The only evidence for a structural change as a 
function of fleld is a reorientation of the triangular lat- 
tice [|| starting at fields of around 3 T. In Bi2212 the 
vortex lattice cannot be observed using SANS at high 
fields, due to the disorder induced by point pinning Jl6[ . 
It is interesting to notice that the early VL-calculations 
based on d-wave gap functions predict the square 

lattice to be oriented along the Cu-O-Cu bonds ({1, 1} di- 
rection in LSCO), which is exactly what our experimen- 
tal data reveal. However, more recent VL-calculations 
p5| , p6[ favor a square lattice tilted by 45° from the Cu- 
O-Cu bonds. Our results provide an important stimulus 
to resolve these theoretical discrepancies. In principle 
other sources of anisotropy, such as those involving the 
presence of dynamical stripes or charge-density waves 
p8| , could lead to the formation of a square vortex lat- 
tice. However, the plausibility of this explanation is as 
yet unclear, since our measurements were performed in 
the overdoped regime of the HTSC where a homogeneous 
Fermi liquid description of the electronic properties of 
these materials is currently believed to be adequate. 

In conclusion we have used SANS to make the first 
clear observation of the vortex lattice in LSCO and ob- 
tained strong evidence for a field-induced crossover from 
a triangular to a square lattice. This may reflect the in- 
creasing importance of the anisotropic vortex cores in this 
d-wave superconductor, or coupling to other sources of 
anisotropy such as those provided by charge/stripe fluc- 
tuations. It remains a challenge to corroborate detailed 
SANS measurements with other microscopic data in or- 
der to explain the origin of this exotic vortex behaviour. 
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FIG. 1. SANS diffraction patterns taken at T=1.5 K, after 
field cooling from 40 K in B=0.1 T. A background taken at 
T=40 K has been subtracted. One of the {1, 1} in-plane axes 
is aligned at 45 degrees to the horizontal axis in the plane 
perpendicular to the field. In a) the c-axis lies along the 
field direction and the field was kept constant during cooling, 
high-resolution mode: Ao — 14.5 A. In b) the c-axis has 
been rotated 10° (about the vertical axis) away from the field 
direction, and while cooling the field was oscillated about its 
mean value, Ao = 8 A. 
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FIG. 2. SANS diffraction patterns obtained by subtract- 
ing the zero-field data at 5 K from the data taken at 5 K for 
a field B=0.8 T applied parallel to the beam, Aq = 8 A. One 
of the {1, 1} in-plane axes is aligned at 32 degrees to the hori- 
zontal axis in the plane perpendicular to the field. The c-axis 
of the sample is either a) parallel to the field or b) rotated 
30 degrees away from it. The result of the two-dimensional 
fits of the diffraction patterns are shown in a) by the white 
square and in b) by the white rhomboid which illustrates the 
distortion induced by the rotation of the c-axis away from 
direction of the applied magnetic field. 
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FIG. 3. a) Tangential average of the neutron signal for 
B=0.2 T (circles), 0.5 T (squares) and 0.8 T (triangles) at 
T=5 K. The intensities have been rescaled to allow a com- 
mon vertical axis. The shift in q illustrates the expected de- 
pendence of the vortex lattice spacing on increasing magnetic 
field. The instrumental resolution (FWHM) is given by the 
horizontal segments, b) The black circles show the field de- 
pendence of cr = {2n/q)^B/$o at T=1.5 K. The horizontal 
lines represent the expected values of a for triangular and 
square lattices. The open circles show the intensity ratio of 
sectors containing the {1,0} and {1, 1} directions. 
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